Mon. Not. R. Astron. Soc. 000, [THS] (2007) Printed 14 March 2009 (MN WT^ stylo file v2.2) 



0^ 

o 

O 



The determination of the rotation period and magnetic 
field geometry of the strongly magnetic roAp star 
HD 154708^ 



S. Hubrigi'2f, G. Mathys^, D. W. Kurtz^, M. Scholler^, V. G. Elkin^, H. F. Henrichs^ 

^European Southern Observatory, Casilla 19001, Santiago, Chile 

^ Astrophysical Institute Potsdam, An der Sternwarte 16, D-14482 Potsdam, Germany 

'^Centre for Astrophysics, University of Central Lancashire, Preston PRl 2HE, UK 

'^European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching hei Miinchen, Germany 

^Astronomical Institute, University of Amsterdam, Kruislaan 403, 1098 S J Amsterdam, Netherlands 



6 



> 

in 

(N 

rn 
O 
On 
O 



Accepted 2007 Enero 99. Received 2007 Encro 98 



ABSTRACT 

We obtained thirteen spectropolarimetric observations of the strongly magnetic 
rapidly oscillating Ap star HD 154708 over three months with the multi-mode instru- 
ment FORSl, installed at the 8-m Kueyen telescope of the VLT. These observations 
have been used for the determination of the rotation period of P = 5.3666 ± 0.0007 d. 
Using stellar fundamental parameters and the longitudinal magnetic field phase curve, 
we briefly discuss the magnetic field geometry. The star is observed nearly pole-on and 
the magnetic field geometry can be described by a centred dipole with a surface po- 
lar magnetic field strength between 26.1 and 28.8 kG and an inclination of the 
magnetic axis to the rotation axis in the range 22.5° to 35.5°. 
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1 INTRODUCTION 

The magnetic A and F stars have global magnetic fields 
with surface field strengths typically of a few kG to tens 
of kG. A subset of them, the rapidly oscillating Ap (roAp) 
stars, pulsate in high radial overtone magnetically distorted 
dipole and quadrupole modes w ith periods in th e range 
5.65 - 21 min (see, e.g.. Table 1 of lKurtz et "al]|2006h . These 
stars thus oflFer the opportunity to observe the interaction 
of p modes with strong magnetic fields as can be done for 
no other star but the sun. The kG-strength magnetic fields 
in sunspots dissipate energy in solar p-modes and shift their 
phases significantly, as is clearly seen in both observational 
and theoretical local helioseismology. While it is not pos- 
sible to perform local asteroseismology in the same detail 
as for the sun - for obvious reasons of angular resolution 
- some unique properties of the roAp stars have the poten- 
tial to allow 3D mapping of their pulsation modes, magnetic 
field geometries and abundance distributions within the ob- 
servable layers of their atmospheres, hence allow a detailed 
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study of the interaction of the pulsations with the magnetic 
fields. 

The roAp stars with very strong magnetic fields are 
of special interest because the efi'ect of the magnetic field 
on the structure of their atmospheres can be studied 
with the greatest detail and accuracy. The cool roAp star 
HD 154708 possesses the strongest magnetic field among the 
roAp stars, wi th a mean magnet ic field modulus (B) — 
24.5 ± 1.0 kG jHubrig et al.l 120051 ). which is about a fac- 
tor of t hree greater than that of HD 166473, (B) ~ 5.5- 
9.0 kG (|Mathvs. Kurtz, fc Elkinll2007^ . the roAp star with 
the second-strongest magnetic field. For all of the roAp stars 
the pulsations observed in radial velocity variations of lines 
of the rare earth elements arise at atmospheric heights where 
the magnetic pressure exceeds the gas pressure and the 
Alfven velocity is greater than the acoustic velocity, so the 
pulsations are primarily magnetic with an acoustic compo- 
nent. There have been many theoretical investigations of the 
roAp stars discussing, e.g., mode conversion and reflection, 
contributions of horizontal and vertical motions to the ob- 
served pulsational radial velocities, an d mode excitation and 
suppression in these stars (see, e.g., ISousa fc Cunhal [20081 
and[Saio 2005,, and references therein). 

All of the theoretical models are for atmospheres with 
magnetic field strengths less than 10 kG. Thus the discovery 



2 S. Hubrig et al. 



of pulsation in HD 154708 in the presence of its 24.5-kG field 
charts new territory. With its field strength more than 15 
times the 1.5 kG typically observed in sunspots, HD 154708 
by far shows the most extreme observab le case of pulsatio n 
in the presence of a strong magnetic field (|Kurtz et al.ll200^ . 
We thus wish to know as much about it as possible, and for 
obliquely rotating and pulsating Ap stars it is imperative 
to determine their rotation periods and magnetic variations 
with rotation to extract important geometrical information. 

The atmospheric parameters of HD 154708, T^b = 
6800 K and logg = 4.11, have been determined us- 
ing the B2— G temperature c alibration of Geneva pho- 
tometry (iHauck fc NorthI 1 19931. Egn. 2) and the Hippar- 
cos parallax |Hubris: et aL |2005|): Stromgr en photometry 
(|Napiwotzki. Schoenberner. fc Wenskell993l ) suggests a hot- 
ter temperature, Teff = 7500 K. The magnetic field in this 
star was discovered during our systematic study of magnetic 
fields in chemically peculiar Ap and Bp stars with the multi- 
mode instrument FORS 1 installed at the 8-m Kueyen tele- 
scope (|Hubrig et al.|[2005l ) . Among the whole class of Ap and 
Bp stars, the only star known at that time to have a stronger 
dipolar magnetic fie ld than HD 154 708 was the much hot- 
ter (Teff f^l4 500K: [Leckronel [l973 ) Bp Si star HD215441 

i Babcock's star ) with a mean field modu l us of (-B) = 34 kG 
Babcock|[l960l ). Now lFrevhammer et al] (120081 ) have found 
a field of (B) = 30.29 ± 0.08 kG in HD 75049, making this 
the second-strongest dipolar field known among the upper 
main sequence chemically peculiar stars. 

The prese nce of pulsat i ons w ith a period of 8min was 
discovered bv iKurtz et al.l (|2006l ) using 2.5 h of high time 
resolution UVES (Ultraviolet and Visual Echelle Spectro- 
graph) spectra. They found that the radial velocity ampli- 
tudes in the rare earth element lines of Nd ll, Nd ill, and Pr ill 
were unusually low, suggesting that roAp stars with stronger 
magnetic fields have lower pulsation amplitudes. Further 
study of the pulsation of HD 154708 over its now-known 
5.3666-d rotation cycle will illuminate this suggestion, as 
would theoretical studies of models with field strengths up 
to 30 kG. 

A recent abundance analysis of HD 154708 was based 
on a high spectral re solution UVES spectrum obtain ed on 
September 18, 2005 (|Nesvacil. Hubrig. fc Khanll2008l ). The 
atmospheric chemical composition of this star was shown to 
be typical of cool Ap stars with a significant ionization dis- 
equilibrium for the first and second rare earth ions, which 
is co mmonly observed in the at mospheres of pulsating roAp 
stars (|Rvabchikova et al.|[20o3 ). The model atmosphere wa s 
calculated with the LLModels code (jShulvak et al.ll2004h . 
Light elements, as well as Ti, Fe, and Ni were found un- 
derabundant whereas Sc and Cr showed solar abundances. 
Heavy elements and rare earth elements were strongly over- 
abundant. Ndlll is 0.86 dex overabundant with respect to 
Ndll and Prill was overabundant by 1.11 dex compared to 
Prll^ 

ICowlev et ai] (|2007D studied the presence of heavy Ca 
isotopes in HD 154708 using the infrared triplet lines of 
Can. They measured isotopic shifts of -1-0.14 A, -1-0.26 A, 
and -1-0.29 A for A8498, 8542, and 8662, respectively. How- 
ever, the maximum isotopic shifts between the isotopes ''"Ca 
and **Ca for the Can triplet lines is of the order of 0.2 A. 
Possibly these differences are caused by other effects than 
isotopic shifts - most likely, blending with highly displaced 



Table 1. Magnetic field measurements of HD 154708 with 
FORSl. All quoted errors are la uncertainties. In the last three 
lines w e list the earlier measurements published in iHubrig et al.l 
1I2OO6I) . 
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54209.3072 


7941±15 


7711±41 


54215.2792 


7557-1-1 5 


7301 +49 


54223.1734 


6765±15 


6174±41 


54238.2840 


6427±23 


6062±65 


54247.1726 
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54254.1219 


6497±28 


6103±80 


54258.2500 


7462±20 


7330±59 


54270.3097 


6479±15 


6169±46 


54279.1667 


7884±15 


7710±44 


54287.2260 


6326±59 


5718±169 


54297.3028 


6450±28 


6211±83 


54305.1536 


7873±18 


7492±49 


54307.0214 


7032±18 


6884±41 


53120.376 




7530±54 


53487.302 




5764±25 


53519.344 




5819±52 



Zeeman components. The magnetic field of HD 154708 is so 
strong that a number of spectral lines are distorted beyond 
recognition by the partial Paschen-Back effect. 

The magnetic field geometry in HD 154708 has not pre- 
viously been studied. In this work we present for the first 
time a mean longitudinal magnetic field measurement series, 
obtained with the multi-mode instrument FORS 1 installed 
at the 8-m Kueyen telescope. The data are used to obtain 
the rotation period and to put constraints on the magnetic 
field geometry. This is the first, necessary step for a more 
detailed study of this remarkable star. 



2 OBSERVATIONS 

The observations were carried out in 2007 April - July in 
service mode. Using the narrowest available slit width of 0'.'4 
the achieved spectral resolving power of the FORS 1 spec- 
tra obtained with the GRISM 600B was about 2000. These 
observations made use of a new mosaic detector with blue 
optimized E2V chips, which was installed in FORS 1 at the 
beginning of April 2007. It has a pixel size of 15 /im (com- 
pared to 24 fim for the previous Tektronix chip) and higher 
efficiency in the wavelength range below 6000 A. With the 
new mosaic detector and the GRISM 600B, we are also able 
to cover a much larger spectral range, from 3250 to 6215 A. 
A detailed description of the assessment of the longitudinal 
magnetic field measureme nts using FORS 1 is presented in 
our previous papers (e.g., iHubrig et aLll2004al lbl. and refer- 
ences therein). We repeat here the major steps of the mag- 
netic field determination. The mean longitudinal magnetic 
field, (-Bz), was derived using 

\2 



V_ 
I 



gaaeX 1 dl 



(1) 



47rmec2 I dX 

where V is the Stokes parameter which measures the circular 
polarization, I is the intensity in the unpolarized spectrum, 
QcS is the effective Lande factor, e is the electron charge, A is 
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the wavelength, me the electron mass, c the speed of light, 
d//dA is the derivative of Stokes /, and (Bz) is the mean 
longitudinal magnetic field. To minimize the cross-talk ef- 
fect, we executed a sequence of spectropolarimetric observa- 
tions at different position angles of the retarder waveplate, 
-1-45—45, -f45— 45, +45—45, etc., and calculated the values 
V/I using: 



V_ 
I 



r - r 
f° + r 



r - r 



(2) 



where a denotes the position angle of the retarder wave- 
plate and /° and axe ordinary and extraordinary beams, 
respectively. Stokes / values were obtained from the sum 
of the ordinary and extraordinary beams. To derive (-Bz), a 
least-squares technique was used to minimize the expression 



{B^)x,-hf 



where, for each spectral point i, yi — {V/I)i 



(3) 



^47rm~~c2" (1/^ ^ d//dA)i, and fe is a constant term that, 
assuming that Eq.[T]is correct, approximates the fraction of 
instrumental polarization not removed after the application 
of Eq. [2]to the observations. In our calculations we assumed 
a Lande factor gctt — 1 for hydrogen lines and QcB = 1.25 
for metal lines. The average Lande factor for metal lines 
is obtained using individual Lande factors for the observed 
sp ectral lines and is close to the value gfoff = 1.23 published 
bv lRomanvukI (|l984l ). Furthermore, in our reduction proce- 
dure the spectral regions containing telluric lines have been 
excluded in the measurements of the magnetic field. The 
errors of the measurements of the polarization have been 
determined from photon statistics and have been converted 
to errors of field measurements. Longitudinal magnetic fields 
were measured in two ways: using only the absorption hy- 
drogen Balmer lines, or using the entire spectrum including 
all available absorption lines. Due to the low resolution of 
the FORSl polarimetric spectra and due to the fact that the 
major contribution to the measured magnetic field strength 
comes from the hydrogen lines, the weak-field approxima- 
tion can be reasonably adopted for magnetic field strengths 
up to a few tens of kG. The results of our magnetic field 
measurements are presented in Table [1] In the first column 
we provide the modified Julian dates at the middle of the 
exposures. The measured mean longitudinal magnetic field 
(-Bz) using the whole spectrum is presented in Col. 2. The 
measured mean longitudinal magnetic field (Bz) using all 
hydrogen lines is listed in Col. 3. All quoted errors are la 
uncertainties. In the la st three lines we list the earlier mea- 
surements published in iHubrig et all (|2006l ). 



3 PERIOD DETERMINATION 

An initial frequency analysis was performed on the longitu- 
dinal magnetic field measurements given in Col. 2 of Table [1] 
using the discrete Fourier transform programme of iKurtj 
(| 19851 ). The resulting amplitude spectrum clearly shows a 
dominant peak with an equivalent period near 5.3 d, but 
the 13 measurements are not uniformly sampled over the 
rotation cycle for this period. We thus used the more ap- 
propriate frequency analysis technique of fitting a sinusoid 



by linear least-squares for a range of frequencies. Fig. [T] 
shows the unambiguous result. The best-fitting frequency 
of 0.18627 d~^ was then fitted by both linear least-squares 
and nonlinear least-squares sinusoidal fitting routines. The 
latter allows an estimate of the frequency uncertainty, and 
the former gives estimates of the amplitude and phase er- 
ror of the fit. The resulting best fit gives the rotational fre- 
quency /rot ~ 0.1863±0.0007d~^, for a rotational period of 
Prot = 5.367±0.020 d. The derived ephemeris for this period 



MJD54257.26 ± 0.03 + 5.367 ± 0.020B 



(4) 



To refine the period determination we combined the data in 
Col. 3 of Table[T]for the hydrogen lines with three previously 
publishe d measurements mad e using the same technique ob- 
tained bv lHubrig et al.l (|2006l ). This data set has a time span 
of 3 y, hence higher frequency resolution than the frequency 
analysis of the data obtained in a single season presented 
in Fig. [T] However, with only a few data points from pre- 
vious years, the alias problem must be examined carefully. 
Because of the sparse distribution of the data, a frequency 
analysis was performed by fitting sinusoids by least-squares 
over an appropriate range of frequencies. As can be seen in 
Fig. [2j there is some alias ambiguity, but it is not severe and 
a best-fitting frequency can be selected. As long as this is 
the correct alias, the rotation period is therefore refined to 
Prot = 5.366 ± 0.002 d for HD 154708. For this determina- 
tion equal weights were used. To increase the accuracy of 
the rotation period even further, we used a non-linear least- 
squares fit of the function (Bz) -I- A(^B^)Cos{2-K{t / P + 4>)) to 
the same 16 data points with weights equal to l/cr'^, with the 
Icr uncertainties as listed in Col. 3 of Table [T] We used the 
value P = 5.366 found as the starting value for the fit and 
find Prot = 5.3666 ± 0.0007 d. We consider this rotation pe- 
riod as the most accurate. A similar fit was performed on the 
13 data points obtained using the entire spectrum including 
all available absorption lines, but we kept in this case the pe- 
riod fixed. The resulting final values are (Bz) — 7199 ± 11, 
= 789 ± 13 and <h = 0.774 ± 0.003 with a reduced 



A 



<Sz> 



X = 3.6. This implies for the derived ephemeris for the 
maximum value of the magnetic field: 



(Bz 



MJD54257.24 ± 0.02 + 5.3666 ± 0.0007B 



(5) 



where the reference date is at the maximum closest to the 
middle of the observations presented in this paper. 

The corresponding phase diagram, including the best 
sinusoidal fit, is shown in Fig. O along with the residuals 
(lower panel). Note that the size of the error bars is com- 
parable to the size of the dots representing the individual 
meeisurements. The scatter of the latter about the best-fit 
curve is also of the same order: this is a strong indication 
that the la uncertainties given in Table[T]are goo d estimates 
of the actual random errors of the measurements (|Landstreetl 
|1980D . To our knowledge, the present data represent the first 
set of mean longitudinal magnetic field measurements of a 
single star obtained with FORS 1 with a sufficiently uniform 
phase coverage to establish this important result. The mea- 
surement with the largest sigma {{Bz)^^^ — 6326 ±59 G) was 
obtained in weather conditions classified as "thick clouds", 
where the guide star was frequently lost, and thus was re- 
peated by the service observer a couple of nights later. With- 
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Figure 1. For the data in Col. 2 of Table [T] this plot shows 
the standard deviation of one measurement with respect to a 
sinusoidal fit for a range of frequencies. There is a clear best fit 
for a frequency of 0.1863d^^, or P = 5.367d. 




0.000 0.125 0.250 0.375 0.500 
frequency 




0.160 0.172 0.185 0.197 0.210 
frequency 

Figure 2. For the data in Col. 3 of Table[T] plus three previously 
published measurements of the magnetic field of HD 154708 for 
the hydrogen lines, these plots show the standard deviation of 
one measurement with respect to a sinusoidal fit for a range of 
frequencies. The top panel shows that the best-fitting frequency 
is close to that determined in Fig.[T] and the bottom panel shows 
at higher resolution that there is a best-fitting frequency with 
some alias ambiguities that fit less well. The best-fitting frequency 
is /rot = 0.18635 ± O.OOOOSd-i, or Prot = 5.367 ± 0.002d. A 
weighted nonlinear least squares fit finally gives Prot = 5.3666 ± 
0.0007 d (see eq. O and discussion in the text. 
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Figure 3. Phase diagram for the magnetic field measurements 
using all lines (Col. 2 of Table [TJ for the most most accurate 
period, Prot = 5.366 d. The variation has a mean of (Pz) = 7199± 



lie and an amplitude of A 



789 ± 13 G. The residuals 



(Observed - Calculated) for the best fit is shown in the lower 
panel. The deviations are mostly of the same order as the error 
bars, and no systematic trends are obvious, which justifies a single 
sinusoid as a fit function. For the point near phase 0.6, see text. 



out taking into account this measurement we obtain a re- 



duced X = 2.6. 



4 DISCUSSION 

A rather short rotation period of 5.3666 days is obtained for 
HD 154708 based on the new thirteen low resolution spec- 
tropolarimetric ob servations. This is fully consistent with 
the result found bv lMathvs et al.l(|l997^ that mean magnetic 
field moduli in excess of 7.5 kG are only found in stars with 
rotation periods shorter than 150 days. While many mag- 
netic Ap stars have very long rotation periods, up to several 
tens of years, or even close to 1 century in some cases, none 
of t hem have extr e mely strong fields such as HD 154708. 

iHubrig et all (|2005h found for HD 154708 a very low 
value for vsini. Our recent fit to the observed spectral line 
profile of the magnetically insensitive line Fe I A 5434.5 re- 
sulted in vsini = 3.5 ± 0.5kms~^ . This is consistent with 
the va lue of 4kms~^ obtained by lElkin. Kurtz, fc Mathvsl 
(|2008l ). 

Knowing the p osition of the star in the H-R diagram 
(|Hubrig et al.ll2005l '). we determine the radius of the star to 
be in the range R = 1.7 — 2.0 Rq. The equatorial rotation 
velocity is given by Uoq ~ 50.6R/P, where R is the stellar 
radius in solar units and P the period in days. From the 
measured v sin i values and Vcq determined using the known 
rotation period and the radius computed from the luminos- 
ity and effective temperature, the inclination of the stellar 
rotation axis of HD 154708 can be determined. We obtain 
'^cq = 16 — 18.9 kms~^, leading to a value for the inclination 
angle in the range i = 9° — 14.5°. 

Our measurements exhibit a smooth, single-wave vari- 
ation of the longitudinal magnetic field during the stellar 
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Figure 4. One spectrum of HD 154708 was obtained on May 20, 
2005 at phase 0.72 and is indicated by the solid line. The second 
spectrum was observed on September 18, 2005 at phase 0.37 and 
is presented by the dotted line. The third spectrum was observed 
on July 23, 2006 at phase 0.76 and is presented by the dashed line. 
No obvious changes arc detectable in the split Zeeman structure 
of the presented spectral lines. 



rotation cycle. These observations can be considered as ev- 
idence for a dominant dipolar contribution to the magnetic 
field topology. Based on the assump tion tha t the m agnetic 
Ap stars are oblique dipole rotators, IPrestonI (| 19671 ) defined 



cos /3 cos i — sin /3 sin i 



{Bz ) '^'^^ cos f3 cos i + sin (3 sin i 
so that the obliquity angle, (3 is given by 

1 - 



P = arctan 



cot I 



(6) 



(7) 



From the values determined above, (Bz)™'"' = 7988 ± 
17 G and (Bz)""'" = 6410 ± 17 G, we find r = 0.802 ± 0.004, 
which for i = 9° — 14.5° leads to a magnetic obliquity 
in the range of (3 — 22.5° — 35.5°. This value is close 
to the low end of th e rang e of obliquities determined by 
iLandstreet fc Mathvd (|2000l ) for magnetic Ap stars with ro- 
tation periods shorter than 25 days. 

For a tilted, centred magnetic dipole and assuming a 
limb darkening parameter of it = 0.5, the surface polar 
field strength is Bd ^ 3.23 (Bz)™'"' (|Prestonlll967l ). i.e. 
Bd ^ 25.8 kG for HD 154708. Using the values determined 
for i and /3, we find Bd=26.1kG for i=9° and /3=35.5° and 
Bd=28.8kG for i=14.5° and /3=22.5°. The measured mag- 
netic field modulus of ^2 4.5 kG using magnetically resolved 
lines l|Hubrig et al.ll2005l ) is of the same order as that esti- 
mated for the presented geometry of the magnetic field. 

The limb-darkening coefficient is a function of wave- 
length. Since no realistic model atmosphere exists yet for the 
roAp stars with extremely strong magnetic fields, the actual 
mean limb-darkening coefficient for the spectral range we ob- 
served for HD 154708 is uncertain. We note that the choice 



of other mean limb-darkening coefficients would extend the 
range of possible Bd values. 

The high spectral resolution UVES spectra of 
HD 154708 at three different dates at a resolving power of 
A/ A A ~ 1.1 X 10''' available in the ESQ archive show almost 
the same appearance of magnetically resolved Zeeman com- 
ponents. In Fig. U we present these three spectra where no 
significant changes in the split line structures is visible, in- 
dicating no noticeable change in the magnetic field modulus 
between these effectively two rotation phases (since two of 
the spectra happen to be at nearly the same rotation phase). 
The spectra obtained on May 20, 2005 at phase 0.72 and on 
July 23, 2006 at phase 0.76 are average spectra of the time 
series used to study rapid oscillations in this star. This does 
not represent a strong constraint: in more than half of the 
known Ap stars with resolved magnetically split lines, the 
ratio between the extrema of the mean field modulus is lower 
than 1.1, and in 75% of them, it does not exceed 1.3. The 
three measurements of {B) that are currently available are 
insufficient to define its variation, or to check its consistency 
with the deriv ed value of the r otation period. 

[Lan dstrcct fc Mathvsl (|200Ql found a predominance of 
small values of /3 for stars with Prot > 25 d. The only star 
with a comparable field strength in their sample is Bab- 
cock's star HD 215441, which has i and /3 angles of the 
order of ~30°, i.e. similar to the i and /3 determined for 
HD 154708. Babcock's star is much hotter and more mas- 
sive than HD 154708, and its rotation period is somewhat 
longer at (P^ot = 9.5 d) than that of HD 154708. It is pre- 
mature to make any statistical assertion about the efficiency 
of magnetic braking and field geometry evolution just from 
the knowled ge of these two strongly ma g netic stars. On the 
other hand. iHubrig. North, fc SchoUeil (|2007l ') studied the 
behaviour of the angle /3 in Ap stars of different masses and 
ages. The obliquity angle distribution as inferred from the 
distribution of r-values appears random at the time mag- 
netic stars become observable on the H-R diagram. After 
only a short time spent on the main sequence, the obliquity 
angle /3 then tends to reach values close to either 90° or 0° 
for M < 3M0. 

As we already mentioned above, due to the presence 
of the very strong magnetic field in the atmosphere of 
HD 154708, its spectrum is very complex and many spec- 
tral lines need a detailed Paschen-Back treatment. On the 
other hand, due to the low value of v sin i and the relatively 
short rotation period, this star is one of the most suitable 
targets to study various atmospheric effects that interact 
with a strong magnetic field. 
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